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Abstract 

The positions of atoms forming a carbon nanotube are usually described by using a 
system of generators of the symmetry group. Each atomic position corresponds to 
an element of the set Z x {0, 1, n} x {0, 1}, where n is a natural number depending 
on the considered nanotube. We obtain an alternate rather different description by 
starting from a description of the honeycomb lattice in terms of Miller indices. In our 
mathematical model which is a factor space defined by an equivalence relation in the 
set {(^0,^1,^2) £ Z 3 I vo + v 1 + v 2 S {0, 1}} the neighbours of an atomic position 
can be described in a simpler way, and the mathematical objects with geometric or 
physical significance have a simpler and more symmetric form. 
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1 Introduction 



The carbon nanotubes, discovered by Iijima in 1991, have several remarkable 
physical properties (geometry-dependent electronic transport from metallic 
to semiconducting with narrow and moderate band gaps, record-high elastic 
modulus, light weight) and many potential applications (molecular electronic 
devices, fiber reinforcement technologies, flat displays, carbon-based nanotips). 
Extensive experimental and theoretical investigations have been carried out 
on the mechanical and electronic properties of these novel fibers. 

The structure of a single-wall carbon nanotube observed by scanning tunnel- 
ing microscopy [14] can be visualized as the structure obtained by rolling a 
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graphene sheet such that the endpoints of a translation vector are folded one 
onto the other. The geometric and physical properties of the obtained carbon 
nanotube depend on this vector, called the chirality of the tubule. The posi- 
tion of the atoms forming a carbon nanotube are usually described by using 
a system of generators of the corresponding symmetry group. Our purpose 
is to present an alternate mathematical model obtained by starting from a 
three-axes description (that is, a description in terms of Miller indices [13]) of 
the honeycomb lattice. 



2 Honeycomb lattice in a three-axes description 



The vectors e = (2/y/E,o), e x = (-1/VE, l/v 7 ^) , e 2 = (-1/^6, -1/v^) 
corresponding to the vertices of an equilateral triangle form a system of co- 
herent vectors [2], that is, for any v = (v x ,v y ), u = (u x ,u y ) G IR 2 we have 

2 2 2 

v = J2( v > e i) e i (v,u) = '%2(v,e i )(u,e i ) \\v\\ 2 = J2(v, e*) 2 (1) 

i=0 i=0 i=0 



where (, ) is the usual scalar product. The numbers vq = (v,eo), v\ = (v,e\), 
v 2 = (v, e 2 ) satisfy the relation v + V\ + v 2 — and can be regarded as the 
canonical coordinates of v with respect to the system of vectors {e ,ei,e 2 }. 
The space 

£ = {(u ,u u u 2 ) E R 3 | M + ui +u 2 = 0} (2) 



becomes in this way a mathematical model for the geometric plane. The corre- 
spondence between this description and the usual one is given by the isometry 

2 

J : R 2 — > £ Xv = (v ,vi,v 2 ) I _1 ( M o, Ui, u 2 ) = Uiei - ( 3 ) 

i=0 



The representation of a vector v G 1R 2 as a linear combination of eo, ei, e 2 is 
not unique. More exactly, we have 

2 

v = ^2viei (v Q , v 1 , v 2 ) G {(v + a, v 1 + a,v 2 + a) \ a G R}. 

i=0 
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Fig. 1. The honeycomb lattice and the partition defined by a vector c € T. 

For each v G M 2 we denote by (vo,vi,v 2 ) (or simply by v) an element of M 3 
such that v = voeo + viei + v 2 e 2 . One can verify by direct computation that 

2 2 2 

(u, v) = UiVi = UiVi = u iVi (4) 

i=0 i=0 i=0 

for any u,v Gl 2 . 

The points of the plane M 2 corresponding to the elements of the set 

C = [v = (v , vi, v 2 ) G Z 3 | v + Vi + v 2 G {0, 1}} (5) 

are distinct and form [3] the honeycomb lattice 



2 



{ i=0 



v = (vo,v 1 ,v 2 ) G £ } . (6) 



The bijection £ — ► A : v t— > J2f=o allows us to describe A by using £. 
The nearest neighbours of v are 

v° = (v + u(v), vi,v 2 ) v 1 = (v , vi + v(v), v 2 ) v 2 = (v ,v 1 ,v 2 + v{v)) 

where v{v) = (— 1) V0+V1 +v ' 2 . The six points = (V)- 7 corresponding to % ^ j 
are the next-to-nearest neighbours of v, and one can remark that v 11 = v, 

d : £ x £ — >N d(v, u) = \vo — uq\ + \v± — ui| + \v 2 — u 2 \ (7) 
is a distance on £, and a point w is a neighbour of order / of f if d(v , u) = I. 
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We have £ = T U (T + ■&), where = (1, 0, 0) and 

T ={v = (v , v u v 2 ) E Z 3 | v + Vl + v 2 = 0} = £ n Z 3 . (8) 

The symmetry group G of the honeycomb lattice coincides with the group of 
all the isometries of the metric space (£, d), and is generated by the transfor- 
mations [4] 

a : £ — > £ cr(v , vi, v 2 ) = (vi, v 2 , v ) 

g:£ — >£ g(v , v±, v 2 ) = (v , v 2 , vi} (9) 
r : £ — >£ r(v ,v 1 ,v 2 ) = -(v ,v 1 ,v 2 ) +<&. 

The group G contains the subgroup of translations (also denoted by T) 
{£ — >■ £: v\-^v + u\ u<ET} 

generated by a 2 rar and ara 2 r 

a 2 Tar(v , v 1: v 2 ) = (v , vi, v 2 ) + (-1, 1, 0) 
ara 2 T(v , Vi, v 2 ) = (v , Vi, v 2 ) + (-1, 0, 1). 

It is known [12] that the C-C bond length in the case of a graphene sheet 
is 1.44 A. If we use the honeycomb lattice £ as a mathematical model for a 
graphene sheet then we have to multiply the position vectors of the points of 
the lattice by the constant a = 1.44v^6/2 in order to get their lengths in A. 

Consider the Hilbert space (l 2 (£), (,}), where 
l 2 (£) = lil) : £ — > C 



£|^)| 2 <oc} 

v&C ) 



vec 



and the unitary representation of G in l 2 {£) defined by 

g : l 2 (£) _> P(£) {gi)){v)=i){g -i v) . (12 ) 
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Let e be a real number, and 70,71,72 be three complex numbers. The linear 
operator 

2 

H:l 2 (C) — >l 2 (£) (Hip)(v) = eip(v) + ^2'y(v,v i )ip(v i ) (13) 

3=0 



where 



7(1* , v 3 ) 



7j if v(v) = 1 
7j if z/(w) = -1 



(14) 



is a self-adjoint operator 

2 



= e^^i( u )^2(v) + £ £7(^)^1 (u)^^) = (fa, Hfa) . 
v&c j=0vec 



The Hamiltonian used in the tight-binding description of 7r bands in 2D 
graphite, with only first-neighbour C — C interaction, has the form (13). 

Theorem 1. For any k = (ko, k\, k 2 ) G £ the real numbers 



E±(k) =e± |7 e ifeoa + 7ie ifcia + 72 e ifc2a l 

belong to the spectrum of H. 
Proof. The function 



ip : C — ► C ip{y) = < 



a for v G T 
(3 for v <ET + $ 



(15) 



where a, /3 are two constants, is invariant under any translation u e T 
</?(i> + w) = 9?(v) for all v <E £. 



The Bloch type function 



(16) 
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Fig. 2. The first Brillouin zone B (left) and the set B c in case n = 4 (right). 

belonging to an extension of the space / 2 (£) satisfies the relation Hipk = Eipk 
if and only if (a, ($) is a solution of the system of equations 

ea + (7 e ife ° a + 7ie ifeia + 7 2 e ifc2a )/3 = Ea 
(7y oe -ifcoa + 7y ie -i*ia + ^ 2 e~^ a ) a + e(3 = E(3. 

This system has non-trivial solutions if and only if 



e — E 7 e ifcoa + 7!e ifcia + 726" 

7 e _ife ° a + 7 1 e _ifcia + 7 2 e _ifc2a e — E 



= 



that is, if and only if E is one of the numbers E±(k). ■ 

The origin on the energy axis is usually chosen such that e = 0. If 70 = 71 = 
72 = 7 is a real positive number then H is a G-invariant self-adjoint operator 
and its spectrum contains for each k G £ the numbers ±E(k), where 

E(k) = 7 | e ifeoa + e ifeia + e ife2a | 

= 3 + 2 cos(A;o — ki)a + 2 cos(A;i — k 2 )a + 2 cos(A;2 — ko)a. (17) 

The relation (17) allows us to extend the function E : £ — > R to a periodic 
function E : M 3 — >R 

2tt 

1 

2vr , \ „ /, , , 2vT 



E(k , fci, k 2 )=E (k + h, k 2 

= E(k ,k 1 + -^,k 2 > j ^E^^h + ^j (18) 



with the property 

E(ko, ki, k 2 ) = E(ko + a, k± + a, k 2 + a) for all a G 
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The corresponding first Brillouin zone is the hexagonal set (figure 2) 



B = {(k ,k 1 ,k 2 )e£ 



3a 3a 



(19) 



(certain points lying on the frontier of B are equivalent). 

The intervals [-37,0] = {—E(k) \ k e B} and [0,37] = {E(k) \ k e B} 
correspond to the valence n and the conduction n* energy bands, respectively. 
The graphene sheet is a conductor since the gap between these bands is null. 

Since E(k) can be written as 



E(k) =7a/ (cos /cqo+cos fcia + cos fc2a) 2 + (sin fc a+sin fcia + sin k2a) 2 



we have E(k) > 0, and the only points of B with E(k) = are 
'2tt 



± 



3a ' 



2tt 
3a 



\ /2tt 2tt\ / 27T 



2tt\ 
3a"/ 



(20) 



that is, the vertices of the Brillouin zone (usually denoted by K) [10]. 

It is known that the Fermi level for a graphene sheet occurs at the K points. 
The function E : R 3 — > R is not differentiable at these points. We have, for 
example, 



lim 



h 2n 



-0 J2-2cos (fcb-g) 
= 7 I™ J 7 2S 



27 lim 



sin 



3a 
fco-(27r/3a) 



whence 



lim 

ko — > 2-7r/3a 
fco > 2n/3a 



E(ko,-&0)-0 



h - 2n 
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lim 

ko — ► 2-7r/3a 
fco < 2n/3a 



(21) 



E(fc ,-^,o)-0 



ft 3^ 



-7- 



The function E : 

(fc) =7 



& is differentiable at any point k with .E'(fc) 7^ 0, and 
-asin(fc — fci)a + asin(fc 2 — fco)a 



dk 



^3+2 cos(fc — fci)a+2 cos(fci — fc 2 )a + 2 cos(fc 2 — fco) a 



(22) 
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etc. The stationary points lying in the Brillouin zone B are 

(0,0,0), ±( 2 A--,--), ±(--X--), ±(--,--,-) 

\3a 3a 3a ) \ 3a 3a 3a) \ 3a 3a 3a) 



that is, the center (a maximum point, denoted by T) and the middle of the 
edges of the Brillouin zone (saddle points, denoted by M). 



3 Chiral single-wall carbon nanotubes 



A single-wall carbon nanotube can be visualised as the structure obtained 
by rolling a graphene sheet such that the endpoints O and A of a translation 

— * 

vector OA are folded one onto the other (figure 1). The geometric and physical 
properties of the obtained nanotube depend on this vector, called the chirality 
of the tubule and represented in our approach by an element c G T . Without 
loss of generality we can restrict us to the tubules with cq > C\ > c 2 . In the 
case C\ = c 2 we have an armchair tubule, and in the case C\ = a zig-zag 
tubule. The nanotubes with ^ C\ ^ c 2 are called chiral nanotubes. The 
diameter of the nanotube of chirality c is ||c||a/7r. 

After the graphene sheet rolling, the points v — 2c, v — c, v, v + c, v + 2c, ... 
are folded one onto the other, for any v — (vq, v±, v 2 ) G C. Thus, each point of 
the set 

[vq, vi, v 2 ] = v + Zc = { (v + jc , v x + jc 1 ,v 2 + jc 2 ) I j EZ} (23) 



describes the same point of the carbon nanotube of chirality c. Each rational 
number is a class of equivalent fractions, called its representatives. In a similar 
way, for each point of a carbon nanotube we have an infinite number of possi- 
bilities to describe it in our model. A mathematical expression containing the 
coordinates of a point is well-defined only if it does not depend on the repre- 
sentative we choose. We describe the atomic positions on a carbon nanotube 
by using the subset 

C c = {[v ,v u v 2 ]eZ 3 /(Zc) | v + v 1 + v 2 G {0;1} }. (24) 

of the factor space 

Z 3 /(Zc) = {[v] = (v ,v u v 2 ) + Zc | v , ui, v 2 G Z } (25) 



as a mathematical model. One can remark that C c is well defined since the 
condition v Q + v i + v 2 G {0; 1} we impose to [v , Vi,v 2 ] does not depend on the 
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representative we choose. Indeed, (vo+jco) + (vi+jci) + (v2+jc2) = f 0+^1+^2, 
for all j G Z. 

Each point [v] G C c has three nearest neighbours, namely, [v°], [v 1 ], [v 2 ], and 
six next-to-nearest neighbours, namely, [v 01 ], [v w ], [v 12 ], [v 21 ], [v 20 ], [v 02 ]. 

A symmetry transformation of the honeycomb lattice C — ► C : v 1— > gv 
defines the symmetry transformation C c — > C c : [v] 1— > [gv] of the carbon 
nanotube C c if 

[ v ] = [ u ] =>. = [gv] 



that is, if 

f — u G Zc =>- gv — gu G Zc. 



Theorem 2. TTie transformations 



9w ■ b c 
r : C c ■ 



9w[v] = [v + w] 
T [v] = [-V + -&] 



(26) 



are symmetry transformation of C c for all w G T. 

Proof. We have v — w G Zc (f + w) — (u + w) = v — u G Zc, and 
v — u G Zc =^ (— f + i?) — ( — u + = u — v E Zc. ■ 

Let n = gcd{co, ci, C2} be the greatest common divisor of c ,ci,C2, and let 
c = nc', that is, c' = co/n, c[ = c\jn and c 2 = C2/n. The transformation 
g c i represents a rotation of angle 2n/n of the nanotube with respect to its 
axis. Since {c\ — C2)co + (02 — c )ci + (c — ci)c 2 = 0, the vector w = (ci — 
C2, C2 — Co, Co — Ci) is orthogonal to c, and the corresponding transformation 
g w is a pure translation, that is, a translation in the direction of the nanotube 
symmetry axis. The vector b = (1/72.) (ci — c 2 , c 2 — c , c — Ci), where 



H = gcd{ci - c 2 , c 2 - c , c - ci} = < 



n if ci - c' 2 ^ 3Z 
3n if c' x — c 2 G 3Z 



(27) 



defines the shortest pure translation of C c . 

From C0 + C1 + C2 = we get (c : - c 2 ) 2 + (c 2 - c ) 2 + (c - c x ) 2 = 3(cg + cf + c 2 ), 
that is, 72 2 ||&|| 2 = 3||c|| 2 , whence 

g= ^(c 2 + c 2 + c 2 ) GnZ. (28) 
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For any w G T the projections of w on c and b can be written as 

(w,c) ( ci-co c 2 -c \c (w,b) ( c 2 cA 6 

w c = r 1 ^ - + " 2 ^H i w r 1 « " W2 n) q' 



where g' = q/n. It is well-known that in the case of two integer numbers 
77, /i G Z with gcd{?7,/i} = 1 there exist 6 Z with ar] + Pfj, = 1. Since 
gcd{(ci — cq)/1Z, (c 2 — Co)/??.} = 1 and gcd{c 2 /n, ci/ra} = 1 it follows that the 
projection of T on c is Zc/g and the projection of T on 6 is 1b/q' . Let uj £ T 
be the shortest vector with 



If £ c is a chiral nanotube then its symmetry group G c is generated by the 
transformations g c >, and r (additional symmetry operations, namely, mirror 
and glide planes occur only in the case of armchair and zig-zag nanotubes). 
More than that, for any [v] G C c there exist s G Z, m G {0, 1, ...,n — 1} and 
p G {0, 1} uniquely determined such that 

H=r p ^[0,0,0]. (30) 



The usual description of the atomic positions of the atoms forming a carbon 
nanotube [15,16] is based on this remark, and the set 

{(s,m,p) I s G Z, me {0, 1, n - 1}, p G {0,1}} (31) 



is used as a mathematical model. 

The subgroup G c of G c generated by g c i and 5^ is a commutative index-two 
subgroup, and g™, — I, g% = g^, where I is the transformation I[v] = [v] for all 
[v] G C c . The irreducible representations of G c are one-dimensional and can 
be described in terms of generators as [5] 



T, 



(K,m) 



(9c' 



T( K , m ){gu) 



-i2izm/n 
-\na/q' 



where 



m G {0, 1, ...,n - 1} 

k G [0, l-Jiq'ja) 



(32) 



The irreducible representations of G c can be obtained from the irreducible 
representations of G c by using the index-two subgroup induction. 

Consider the Hilbert space (/ 2 (£ c ), (, )), where 
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1 2 {C C ) = U:C C 



E l^)| 2 <oo 

(V^HE^W^) (33) 



and the unitary representation of G c in / 2 (£ c ) defined by 

g : Z 2 (£ c ) — Z 2 (£ c ) (#)M = VGr>])- (34) 

If 5 is a real number and 70, 7i, 72 are three complex numbers, then the linear 
operator 

H : l\C c ) — Z 2 (£ c ) (^)H = e^H +^ 7 («y') (35) 

i=o 

with 7(f,f J ) defined by (14), is a self-adjoint operator. 
Theorem 3. For any k = (k , ki, k 2 ) G £ satisfying the relation 

(k, c) = k c + fciCi + k 2 c 2 G (27r/a)Z (36) 
£/ie reaZ numbers 

E±{k) =e± | 7o e ifcoa + 7l e ifcia + 72 e ifc2a | (37) 

belong to the spectrum of H. 

Proof. For k G £ satisfying (36) we have 

v-ueZc =^ (p(v) e i{k ' v)a = <p(u)e i{k ' u)a 

and hence the Bloch type function 

ij k :C c ^C ^ k [v]=<p(v)^ a (38) 

is well-defined. If E G {E + (k), E_(k)} then there exists a non-null function of 
this form satisfying the relation Hifj k = Eifj k . ■ 

If To, Ti) 72 are re& l then H is G c -invariant. Indeed, since g^v* = (g^vY , 
g^f- 7 = (g u vy and rt^ = (rt> ) J we have gHg^ 1 = H, for any g G G c . Denoting 
(/c, c) = 2mn/a and (A;, a;) = n/q' we get 

(y c ,-0 fe )H = <p(v) e i{k ' v - c ' )a = e- i{k ' c ' )a ^ k [v] = e- i2mw/n ip k [v] 
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(39) 



that is, the subspace generated by ip k and ip- k is G c - invariant. More than that, 
these relations allow us [5] to classify the eigenstates ip k by using the quantum 
numbers m and k. The set (figure 2) 



keS 



(k,c) = 2mir/a with m G {0, 1, n — 1} 
(k,u) = n/q' with k G [0, 2q'ir/a) 



> . 



(40) 



contains a /c corresponding to each class. The conduction m-band 
{E + (k) | (k,c) =2mir/a, (k,u) e [0,2irq'/a)} 



(41) 



and the valence m-band 

{E_(k) | (k,c) =2mir/a, (k,u) e [0,2nq'/a)} 



(42) 



can easily be determined by using (37). A graphic representation of these 
conduction bands can be found in [5]. 

The relation (36) defines a family of equidistant straight lines orthogonal to 
c with the distance between neighbouring lines equal to 5 = 27r/(a||c||). Since 
the length of the projection of the vector (2ir/ (3a), — 27r/ (3a), 0) on c is 2h(cq — 
Ci)/(3a||c||) = (c — ci)S/3, the K points belong to the straight lines (36) if 
and only if c — C\ G 3Z. 

The Hamiltonian used in the tight-binding description of tt bands in C c , with 
only first-neighbour C — C interaction, has the form (35). Except for very small 
diameter nanotubes, the constants 70, 71, 72 keep almost the same values as 
in the case of a graphene sheet [1,16]. Neglecting the effects of the curvature 
of the graphite sheet, we can assume e = 0, 70 = 71 = 72 = 7 G (0, 00), and 
we get the energy levels 



E±(k) = ±^J3 + 2 cos(A;o — hi) a + 2 cos(fci — k 2 )a + 2cos(k 2 — k )a. (43) 



From the form of the surface E(k) [12], it follows that the gap between the 
valence an conduction bands is given by 

AE C = mm E + (k) - max £_(£;) = 2 mm E(k) (44) 

keB c keB c k£B c 
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and we have AE C = if and only if cq — c\ G 3Z. Therefore [7,16], the nanotube 
C c is a conductor if cq — C\ G 3Z, and a semiconductor if Co — C\ G" 3Z. The 
minimum minfc e g c -£"(&) is achieved [12] at a point lying on a straight line (36) 
in the vicinity of a point K. 

In the case of a magnetic field parallel to the nanotube axis the Hamiltonian 
also has the form (35), but 70 = 7e 1/3c ° a , 71 = 7e 1/3cia , 72 = 7e 1/3c2a , where (3 is 
a real constant describing the magnetic field strength threading the nanotube 
[9]. Choosing e = 0, 7 G (0, 00), we get the energy levels 

E± a9 (k) = ± 7 | e ^o+/3c )a + e i(ki+/fci)a + e i(fc 2 +/3c 2 )a | = + ^ ^ 

for any A; satisfying (36). From this relation it follows that the presence of a 
magnetic field parallel to the nanotube axis can modify drastically the elec- 
tronic properties of the nanotube, and this modification depends essentially 
on the nanotube chirality [9]. 



4 Concluding remarks 

The alternate description presented in this paper offers certain formal advan- 
tages and is significantly different from the usual one. We think that it can 
be used as a complementary description, and may stimulate the interest of 
mathematicians in the fascinating geometry of carbon nanotubes. We have 
re-obtained some known results [5,6,8,9,11,12,16] concerning the chiral carbon 
nanotubes in order to illustrate the proposed approach. 

Since our model is a factor space, we have to verify the independence of the 
representative (t>o, v±, V2) we choose for [vq, i>i, 1*2] in the case of any mathemat- 
ical object we consider on C c . This is not an inconvenience for our approach 
but a useful criterion when we look for mathematical objects with possible 
geometric or physical significance . 
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